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1. Introduction 

We are interested in coalescence which is a widespread phenomenon: it arises in physics, chem- 
istry, astrophysics, biology and mathematics. 

We consider a possibly infinite system of particles, each particle being fully identified by its 
mass ranging in the set of positive real numbers. The only mechanism taken into account is the 
coalescence of two particles with masses x and y into a single one with mass x + y at some given 
rate (the "coagulation kernel") K(x,y) = K(y,x) > 0. 

— We can consider a system of microscopic particles and the following system of differential 
equations for the concentrations Im{x) of particles of mass x = 1, 2, 3, ... at time t £ [0, +00): 

^ x — l +00 

(1-1) dtfi t ( x ) = -j^K{y,x-y)vt{y)vt{x-y) - ix t {x)^K{x,y)ix t {y)- 

The first sum in (1.1) on the right corresponds to coagulation of smaller particles to produce 
one of mass x, whereas the second sum corresponds to removal of particles of mass x as they 
in turn coagulate to produce larger particles. 

Analogous integro-differential equations allow us to consider a continuum of masses x. In 
this case the system can also be described by the concentration /x t (a;) of particles of mass 
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x G (0, +00) at time t G [0, +00). Then jJLt{x) solves a nonlinear equation: 



(1.2) 



d t fj, t (x) = - K(y,x-y)fH(y)fi t (x-y)dy- fi t (x) K(x,y)fj, t (y)dy. 





Equation (1.2) is known as the continuous Smoluchowski coagulation equation and (1.1) is 
its discrete version. 

— When the particles are macroscopic and when the rate of coagulation is not infinitesimal, 
the frame of study of the dynamics of such a system is stochastic. When the initial state 
consists of a finite number of macroscopic particles, the stochastic coalescent obviously exists 
(see [1]) and it is known as the Marcus-Lushnikov process. 

In preceding works several results have been obtained on the existence and uniqueness of weak 
solutions to Smoluchowski's coagulation equation. The general framework was formulated in [15] 
who obtained some remarkable well-poscdncss results. In [7], homogeneous-like kernels are con- 
sidered and it has been seen that the well-poscdncss holds in the class of measures having a finite 
moment of order the degree of homogeneity of the coagulation kernel. 

Aldous [1] presents the Marcus-Lushnikov process as an approximation for the solution of Smolu- 
chowski's equation (see [14, 13] for further information). Since then some results on convergence 
have been obtained in [15] and [10], see also [6]. A class of stochastic algorithms in which the 
number of particles remains constant in time was introduced in [3] and has been extended to the 
discrete coagulation- fragmentation case in [11]. 

We investigate the rate of convergence of the Marcus-Lushnikov process to the solution of the 
Smoluchowski coagulation equation as the number of particles tends to infinity. This problem is 
interesting because on the one hand it has a physical meaning: the Smoluchowski equation is often 
derived by passing to the limit in the Marcus-Lushnikov process, and on the other hand from a 
numerical point of view: this stochastic process can be simulated exactly. Thus it seems natural 
to use it in order to approximate the solution to Smoluchowski's coagulation equation. 

Our study is based on the use of a specific Wasserstein-type distance d\ between the solution 
to Smoluchowski's equation and its stochastic approximation. This distance depends on the ho- 
mogeneity parameter A of the coagulation kernel. This specific distance has been introduced in 
[7] to prove some results on the well-posedness of the Smoluchowski coagulation equation and in 
[5, 8] to study the stochastic coalescent. The result of the present work applies to a family of 
homogeneous-like coagulation kernels. These kernels are of particular importance in applications 
see Table 1 in [1] or the list provided in [7]. 

We point out that since we are using a finite particle system to approximate the evolution in 
time of the solution to the Smoluchowski equation which describes an infinite particle system, it 
is necessary to dispose of a mechanism to construct an initial condition for the Marcus-Lushnikov 
process from a general measure-valued initial condition of Smoluchowski's equation. This initial 
condition needs to satisfy, on the one hand, a convergence condition to assure the convergence of 
the stochastic process to the solution to Smoluchowski's equation for all time t as the number of 
particles grows (the usual condition of weak convergence is replaced by convergence in the sense of 
the distance we use) , and on the other hand it must obey a rate of convergence in order to control 
the overall rate of convergence of such an approximation. 



Very roughly, we consider a homogeneous-like coagulation kernel with degree of homogeneity 
A G (—00, 1] \ {0} (including K(x,y) = (x + y) x ). For (nt)t>o the solution to the corresponding 
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Smoluchowski's equation and for (^™) t >o the corresponding Marcus-Lushnikov process, we prove 
that 

sup E[d x (jx2,fit)} < -p, 

te[0,T] V n 

as soon as no satisfies some technical conditions and for a good choice of the initial state of the 
Marcus-Lushnikov process of the form i SfcLi ^ k ■ We can make the following remarks. 

1. Recalling the Central Limit Theorem (CLT), this rate of convergence seems to be optimal, 
since the convergence of /i™ to /it is a generalized Law of Large Numbers. 

2. In [7] it has been seen that only one moment is demanded to show the well-posedness for the 
Smoluchowski equation. In the present work, we need to demand more moments, but we believe 
that it is very difficult to avoid such conditions. 

3. The only works giving an explicit result on the rate of convergence of the Marcus-Lushnikov 
process toward the solution to Smoluchowski's coagulation equation, known by us, are: 

- Norris [15], who gives an estimate using a "Large Deviations" approach for the discrete case 
(supp(^o) C N). 

- Deaconu, Fournier and Tanre [2] , where a CLT- type result is shown for the discrete case and 
for a bounded coagulation kernel K, furthermore in this work a different particle system is 
used. 

- Kolokoltsov [12], who uses analytic methods of the theory of semigroups applied to the 
Markov infinitesimal generator. He also uses a different distance to ours, namely the author 
uses the topology of the dual to the weighted spaces of continuously differcntiable functions 
or certain weighted Sobolev spaces. He then gives a CLT result for the discrete case with 
a coagulation kernel satisfying K{x 1 y) < c(l + y/x)(l + ^/y) and for the continuous case 
when K is two times differentiable with all its derivatives bounded. Unfortunately the case 
K(x, y) = (x + y) x is excluded for any value of A e (-co, 1] \ {0}. 

Our work thus gives the first result on the rate of convergence covering the continuous case for 
some homogeneous kernels. 

For the case A < we follow the ideas found in [7], but for the case A <E (0, 1] the proof is 
much more difficult and the calculations are faced in a completely different way. Namely we use 
the Ito formula for an approximation of the absolute value function and handle very delicately the 
resulting terms. 

The paper is organized as follows: in Section 2 we give the notation and definitions we use in 
this document, in Section 3 we state our main result. The proof is developed in Sections 4, 5 and 
6. We give also a method to construct an initial condition for the Marcus-Lushnikov process in 
Section 7 and we conclude the document giving some technical details which arc useful all along 
the paper in Appendix A. 

2. Notation, Assumptions and Definitions 

In this section we present our assumptions, give the definition of weak solutions to Smolu- 
chowski's coagulation equation and then we recall the dynamics of the Marcus-Lushnikov process. 

Notation 2.1. — We denote by M + the space of non-negative Radon measures on (0,+oo). For 
a measure /i and a function 4>, we set (fj,(dx) , <f>(x)) = J^°° 4>(x) fj,(dx) . We also define the operator 
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A for all measurable functions <j) : (0, +00) — > R, by 

(2.1) (M)(x, y) = 0(x + 2/) - </>(x) - 0(y) V (x, j/) 6 (0, +00) 2 . 

Finally, we will use the notation x A y = min{a;, y} and x V y — max{x, y} for (x, y) G (0, +00) 2 . 

We consider a coagulation kernel K : (0, +00) x (0, +00) — > [0, +00), symmetric i.e. if(x, y) = 
A'(y, x) for (x, y) G (0, +00) 2 . We further assume it belongs to W 1 ' 00 ((e, 1/e) 2 ) for every e G (0, 1) 
and one of the following conditions V (x, y) G (0, +00) 2 : 

(2.2) A G (-00, 0), A"(.T,y) < K (x + y) x and (.t a + y A ) \d x K{x, y)\ < n 1 x x - 1 y x , 

(2.3) A G (0, 1], A>, y) < k (s + y) x and (.t a A y A ) \8 x K(x, y)\ < KiX A " V, 

(2.4) AG (0,1], K(x,y) < K (xAy) x and (x x Ay x ) \d x K(x,y)\ < K l x x ~ 1 y x , 

for some positive constants kq and k\. We refer to [7] for a list of physical kernels satisfying 
conditions (2.2) and (2.3). Remark that for any A G (— 00, 1] \ {0}, K{x,y) = (x + y) x satisfies 
(2.2) or (2.3). 

Definition 2.2. — Consider A G (-00, 1] \ {0}. For fx G M + , we set: 

(2.5) M x (li) = \ x x n(dx) and M\ = {v G : M\ (y) < +00}. 

Jo 

For \i G A / J + , we set, for x G (0, +00): 

(2.6) F»{x) = l (x>+oo) (y) »(dy) and G"(x) = / l(o,x](lO /*(<*»)■ 

Jo Jo 

We define the distance on Ait as 

r+00 

(2.7) d\(fi,,fi) = / ^"^(a^dx, 

Jo 

w/iere £(x) = C(x) - G^(x) if A G (-00, 0) and £(x) = ^(x) - i^(x) i/AS (0,1]. 

We remark that d> is well-defined on M^. Indeed we have d\(fj,, fi) < tjjMx {p. + fi) for A G 
(—00, 1] \ {0}. See [4] for a deeper study of this distance in the discrete and continuous cases. 

We excluded the case A = for two reasons. First, do is not well-defined on Mq . Next, when 
trying to extend our study to this case, we are not able to obtain a better result than those of 
Kolokoltsov [12]. 

Definition 2.3. — For A G (— 00, 1] \ {0} we introduce the spaces of test functions needed to 
define weak solutions: 

if A G (—00, 0) : H\ = {0 : (0, +00) — s- R such that sup 2 , >0 x~ x \<j)(x)\ < +00} , 

if A G (0, 1] : U\ = {0 : (0, +00) R such that sup 3 . > o(l + x) x \tf>(x)\ < +00} , 

if A G (0,1] : U% = {(/) : (0,+oo) -> R such that sup x>0 x~ x \<j)(x)\ < +00} . 

It is necessary to introduce the space % e x to study the case (2.4). 
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2.1. The Smoluchowki coagulation equation. — The weak formulation of the Smoluchowski 
coagulation equation is given by 

(2.8) j t {m{dx),(l>{x)) = ^(jH{dx)iH{dy),{M)^y)K{x,y)), 

see Notation 2.1. This is a general formulation and it embraces the two previous equations : if 
Ho is discrete (i.e. supp(// ) C N), then this corresponds to the "discrete coagulation equation" 
(1.1), while when fj,o is continuous (i.e. Ho(dx) = [io(x)dx), this corresponds to the "continuous 
coagulation equation" (1.2). Formulation (2.8) is standard, sec [15]. 

Definition 2.4- — Let A G (— oo, 1] \ {0}, a coagulation kernel K satisfying either (2.2), (2.3) 
or (2.4), and fj, m 6 M\. We will then say that (Ht)t>o C Ai + is a (fj, in , K, X)-weak solution to 
Smoluchowski 's equation if the following conditions are verified: 

(i) Ho = H m , 

(ii) the application t 1 — > (fit(dx) , <p(x)) is differentiable on [0,+oo) and satisfies (2.8) for each 

cj> e H x (cases (2.2) and (2.3)) or for each cj> e U\ (case (2.4)), 
(Hi) for all T £ [0, +00) 

(2.9) sup M a (n s ) < +00, 

se[o,T] 

for a = A (cases (2.2) and (2.4)) or for a — 0, 2A (case (2.3)). 

We demand more finite moments of fio than in [7] to assure the convergence of the Marcus- 
Lushnikov process. According to the hypothesis on the kernel (2.2), (2.3) or (2.4) together with 
(2.9) and Lemma A.l, the integrals in the weak formulation (2.8) are absolutely convergent and 
bounded with respect to t € [0, T] for every T. 

Under (2.2) or (2.4), the existence and uniqueness of such weak solutions have been established 
in [7] for any [i m € M\. Under (2.3), the existence and uniqueness of weak solutions satisfying 
(2.9) with a — A have also been checked in [7] for any fi m € Mt. Using furthermore Proposition 
A. 4, we immediately deduce the existence and uniqueness of weak solutions under (2.3), in the 
sense of Definition 2.4, for any /x m <S Mq fl M.^- 

2.2. The Marcus-Lushnikov process. — The Marcus-Lushnikov process describes the stochas- 
tic Markov evolution of a finite particle system of coalescing particles. We consider a coagu- 
lation kernel K and a finite particle system initially consisting of N > 2 particles of masses 
Xi, ■ ■ ■ , xm € (0, +00). We assume that the system evolves according to the following dynamics: 
each pair of particles (of masses x and y) coalesce (i.e. disappears and forms a new particle of 
mass x + y) with a rate proportional to K(x, y). 

Let n G N and we assign to all particles the weight 1/n. We define now rigorously the Marcus- 
Lushnikov process to be used. 

Definition 2.5. — We consider a coagulation kernel K, n 6 N and an initial state /ip = 
~ J2iLi $xi, with xi, ■ ■ ■ , xisr € (0, +00). 

The Marcus-Lushnikov process (/x™)t>o associated with (n, K, fj,o) is a Markov M. + -valued cadlag 
process satisfying: 

(i) (/U")t>o takes its values in j — X)i=i ^m'i ^ — y% > 1. 
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(it) Its infinitesimal generator is given, for all mesurable functions ^ : A4 + — > K and all states 

l<i<j<k 

This process is known to be well-defined and unique, see [1, 15]. We will use the follow- 
ing classical representation of the Marcus-Lushnikov process (see e.g. [5, 8]): there is a Pois- 
son measure J(dt,d(i,j),dz) on [0,+oo) x {(i 7 j) £ N 2 ,i < j} x [0, +oo) with intensity measure 
dt [X)fc<; $(k,i) (d(i,j))] dz, such that for any measurable function (f> : (0, +oo) — > R 

pi p p ~\~ oo -i 

(tf(dx),4>(x)) = (^(dx),4>(x)) + / / / -Ufx i a _+xi_)-tf>(x i B _)-4,(xi_)' 

JO Ji<j Jo 71 L V ' 

(2.10) lj ^ K ( xl _ xi _^t {j <: N{s „ )} J(ds,d{i,j),dz), 

where /i" = ^ EfcH* $x k > N(t) being the (non-increasing) number of particles at time t. 
This can be written using the compensated Poisson measure related to J: 

(tf(dx)A(x)) = (^(dx),<t>(x)) + I f t (^(dx)^(dy),(A ( t>)(x,y)K(x,y))d S 



(2.H) -— / (^(dx),(A^)(x,x)K(x,x))ds 



J(ds,d(i,j),dz), 

where the operator A is defined in (2.1). The third term on the right-hand side is issued from the 
impossibility of coalescence of a particle with itself. 



3. Results 

We state in this section our main result. We also state as a proposition the construction of 
a sequence of initial conditions for the Marcus-Lushnikov processes and finally comment on our 
results. 



Theorem 3.1. — We consider A £ (— oo, 1] \ {0} and a coagulation kernel K satisfying either 
(2.2), (2.3) or (2.4)- Let )iq £ M. + and (/x*)t>o the (no, K, \)-weak solution to Smoluchowski's 
equation. Let /1q be deterministic and of the form ^ E^Li an d denote by (/Lt™)oo the associated 
(n, K, fig) -Marcus-Lushnikov process. Let e > 0. 

• Assume (2.2) or (2.4) and that no belongs to A4~t( n M.2\+'> where e = sgn(X) x e. Then for 
any T > 0, 



E 



sup dx(Ht^t) 
te[o,T] 



< 



d\(fj%,fM)) 



(1 + T)C Xtl 



M x {^)+M 2X+i {^) 



xexp[TC A , £ M A (^ + Ai )], 
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where C\_ £ is a positive constant depending only on X, e and kq, and K\. 
• Assume (2.3) and that ^o € M$ fl M.~ +E where 7 = max{2A, 4A — 1}. Then for any T > 0, 



sup E[d A (A#,^)] < 
te[o,T] 



(1 + r)c v 



l + [M (^ + Mo)] 



+ [M 7+E (/# + /*,,)] 



xexp[TC A , e M A (^ + M o)], 



where C\_ e is a positive constant depending only on X, e, kq and k,\. 

Now we present the proposition giving a d A -approximation of the initial condition. 

Proposition 3.2. — Let X 6 (— 00, 1] \ {0} ; n e N and /xo a non negative Radon measure on 
(0, +00) such that fiQ G -M^ fl M^x- The measure fio is supposed to be either atomless or discrete 
(Wpp(/io) C N). Then, there exists a positive measure fiQ of the form i such that: 

d\(fJ>Q, Mo) < -k, 

where the constant C\ depends only on X and M2a(mo)- We also have 

MM) < M a (fi ), 

for all a < if X 6 (—00, 0) and for all a > i/Ae (0, 1] . Furthermore, if Mo(/io) < +00, i/ien 

JV„ < nMo(Mo). 

The estimate of the parameter iV„ (initial number of particles) may be useful to study the 
numerical cost of the simulation. 

Gathering Theorem 3.1 and Proposition 3.2, we deduce the following statement. 

Corollary 3.3. — We consider X 6 (—00, 1] \ {0} ; e > and a coagulation kernel K satisfying 
either (2.2), (2.3) or (2.4-). Let /io € A4 + be either atomless or discrete (supp([io) C N), and 
(Ht)tG[o,+oo) the (/Iq, K, X)-weak solution to Smoluchowski's equation. Then it is possible to build 

a family of initial conditions Mo = ^ Y^,k=i ^4 such that, for (/i")t>o the corresponding (n, K, /j,q)- 
Marcus-Lushnikov process, 

• under (2.2) or (2.4), if Mo belongs to H M-2\+s> where e = sgn(X) x e, then for any 



T > 0, 



sup d\(nt,nt) 

t6[0,T] 



< 



C T 



where Ct is a positive constant depending only on T , X, e, kq, hi and /io; 
• under (2.3), if Mo e Mq PI M^ +e where 7 = max{2A, 4A — 1}, then for any T > 0, 

sup E[d x (i$,iit)] < 
te[o,T] V n 



where Ct is a positive constant depending only on T, X, e, kq, k\ and mo- 
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This last statement is quite satisfying since it provides a rate of convergence in and it 
applies to a large class of homogeneous kernels presenting singularities for small or large masses. 
We probably require more finite moments than really needed but this docs not seem to be a real 
problem for applications. 

We have followed the ideas found in [7] to prove the case (2.2) and the special case (2.4) of 
Theorem 3.1. The case (2.3) is much more subtle and difficult. For this case we have applied the 
Ito formula and manipulated very carefully each term. By the moment it is not possible to put 
the "sup" into the expectation since it is very important to use the sign of the terms and to take 
advantage of some cancelations. 

Proposition 3.2 presents the proof of the existence of a dA-approximation of a general non- 
negative measure /j,q (we consider measures which are interesting for the Smoluchowski's equa- 
tion) by a discrete measure /Xq (a finite sum of Dirac's deltas) as a construction procedure. This 
construction is very useful from a numerical point of view since it gives a measure that will be set 
as the initial state for the Marcus-Lushnikov process. 



4. Negative Case 

In the whole section, we assume that K satisfies (2.2) for some fixed A E (-co, 0). We fix e > 0, 
and we assume that /zo G M\ H -Mtx—e- We denote by (nt)t>o the unique (/j,q, K, A)-weak solution 
to the Smoluchowski equation. We also consider the (n, K, fig )-Marcus Lushnikov process, for 
some given initial condition /Ltg = i YliLi <W 

We introduce, for t > 0, the quantity E n (t,x) = G^" (x) — G' l4 (x) as defined in (2.6). Wc take 
the test function <f)(v) = lm x ](v). Since sup t)>0 v~ x \4>(v)\ = x~ x < +oo, we deduce that <j> G 7~L\. 
Computing the difference between equations (2.11) and (2.8), we get 

1 /"* 

E n (t,x) = E n (0,x) + -J^ (f^(dv)^(dy)-n 3 (dv)fx s (dy),(Al (0!x] )(v,y)K(v,y))ds 

(4.1) - A- J (fi(dv), (Al ( o, x] ) (v, v)K(v, v))da 

r+°° i / . N 

J n \ X '-^ X '-) 1 j. -..fx x i| ' b v - i 

J(ds, d(i,j),dz). 

We take the absolute value and integrate against x x_1 dx on (0, +oo): 




(4.2) 



dx(j#,iH) < d x {^,no) + A 1 (t) + A 2 {t) + A 3 (t), 
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where 
A 1 (t) 

Mt) 



1 



+00 



-A-1 



2 Jo 

_L /" +0 ° 

+00 

~A-1 



(^(dv)Hs(dy) - n s (dv)n s (dy), (Al (0iX] ) (v,y)K(v,y))ds 



dx, 



.A-l 



(41(0,*]) 



da;, 



+00 



/0 Ji<j "'O 

Now we are going to search for a good upper bound for each term. 
Term A x {t). 

Similarly to [7, Lemma 3.5]. However, in this case we have to argue a little more, since t 
G M "(x) is not (even weakly) differcntiablc due to the jumps of /x". 

The term A\(t), according to the symmetry of the kernel, can be written as: 



M(t) 



2 Jo 



„A-1 



n+00 />+oo 
Jo 



K{v,y) [l^ x] (v + y) - l {0 ,x](v) ~ ^(o,x](y)] 



(4.3) 



(/i™ - fx s ) (dv) (/£ + n s ) (dy)ds 



dx. 



We use the Fubini theorem and Lemma A. 2: 

r-t /> + oo /> + oo 



JO JO 



K(v, y) [1 {0<X ] (v + y)- l(o,»](«) ~ l(o,x](?/)] " Ms) {dv) (/£ + fi s ) (dy)ds 



t r-\-oo r-\-oo 



'0 JO Jo 

r+00 



K(x-y, y)t(o, x ] (v + y) - K(x, y)t( ,x] (v) 



d x K(z,y) [t( ,x]{z +y)~ ^(o,x](z) - ^(o,x](y)] dzj (p" - (M s ) (dv) 

« + A*a) (dy)ds 



t x >y t( 0tX -y)(v) ((J% - Ms) (dv) 



+00 



n+00 
K(x-y,y) 

n+00 
K(x,y) 

d x K(z,y) [l {0 ,x](z + y)~ l(o,z](*0 ~ ^(o,x\(v)\ 



(M? + Ms) (dy)ds 



^(o,x](v) (m™ - Ms) (dv) 



(fj% + fi s ) (dy)ds 



Jo JO 



+00 



l(o,*] («) (M?-Ms) (dv) 



dz (n n s + n s ) (dy)ds 
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t p+oo 

K{x ~ V, y) \^x> y E n (s, x - y)} (/i™ + n s ) (dy)ds 



JO 

K(x, y) [E n (s, x)] (p™ + fi s ) (dy)ds 

- / / / d x K(z,y) [l( 0)X ](z + y)-l(o, x ](z)-l(o iX ](y)] 
Jo Jo Jo 

[E n (s, z)] dz (X 1 + fx s ) (dy)ds. 

According to the bound 

(4-4) \l {0jX] (z + y) -t (0tX] (z) - 1(0,^(2/)! < 2 l(o,x]0* A y), 

and using (2.2), we deduce: 

i.+oo ,.+00 

A^t) < / / / x A -V|£„(s,x-y)|dx(/^+/^)(dy)ds 

^ JO JO Jj/ 
„ z-t /-+oo (.+00 
K ° / 



ir*" 1 ^ + y) x \E n (s, x)\ dx (/£ + fi s ) (dy) ds 



2 Jo Jo Jo 

ft p-\-oo />+oo 

dz (/i™ + // s ) (dy) ds. 







+00 

a ;A_1 l(o, a; ]( z A y)dx 



\d x K(z,y)\\E n (s,z)\ 

'0 Jo Jo 

For the first integral we use the change of variable x H> w+y and (w+y) x_1 (u>+y) A < w x ~ 1 y x . For 
the second integral (x + y) x < y x . Finally for the third integral, we observe that J £ A_1 1(o.:e] (zA 
?/)da; = ^tjt < 2 |^ . Using (2.2) again, this implies 

pt r+00 r+00 

< ds w^E^w^dw y x (n't + us) (dy) 

z Jo Jo Jo 

pt p+CO Z'+OO 

+ _£ / ds / x x -^\E n ( Sl x)\dx y x (fi s ' + n.) (dy) 

* Jo Jo Jo 

pt p+oo p+oo 

+ W\ Jo ^ Jo ^ 1|S " (S ' z)l dz I yX (M " + Ms) 
The resulting bound for A\(t) is: 

(4.5) A x (t)< ( K o + j^)^ d A (A*?,A».)M A (A*?+A*.)da. 



Term A 2 (t). 

We use |(A1 ((M ) (v,v)\ = \l (0iX] (2v) - 21 (0iX] (t;)| = l {0 <t,<f} +21 {|<1 ,<. t} < 2t {v < x} . This 
gives 

^ r+00 pt p+00 

A 2 (t) < - / z*- 1 / / #(u,u)l { „< x} /^(cfo)ds<fo; 
n Jo Jo Jo 

< - / / K (2v) x —^(dv)ds 
n Jo Jo \M 

(4-6) < ^ f MM)ds. 

n \M Jo 

We used (2.2). 
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Term A 3 (t). 



We will bound the expectation of this term using its bracket, for this we consider: 
/ 

i 

\ 



E 



fO Ji<jJQ 
rt l 



(At (0 . x] ) \Xl_,X J s _j 1, ^ K ( xt _ xi _^l {j < N(s _ )} J(ds,d{i,j),dz) 



= E 



4^ 
< — E 
n 



2^ 
< -E 
n 



< — -E 



K ( X* 

E 1 : S) [1(0,,] (Xl+Xl) l (0iX] (XI) 1,0.,] (X$)Yds 



n ^ * — , n 

i<j<N(s) 



E 



KiXIXl) 



u i<j<N{s) 



1(0,,] (XlAXl)ds 



IV 



(fi(dv)fi(dv),K(v,y) [l(o,x](«) + l(o,x](tf)]) ds 
(tf(dv)tf(dy), (v + y) x t (0tX] (v)) ds 



We have used (4.4), a symmetry argument then the bound t(o. x ]{v V y) < l(n tX ](v) + 1( . x ](y) and 
finally (2.2). We consider now the submartingale (absolute value of a martingale): 



S t (x) 



«/i<7 JQ 



(-41(0,,]) (^->^f-J !f < K(x*_,xi_) > [ 1 {j<N(s-)} J(ds,d(i,j),dz) 



According to the Cauchy-Schwartz and Doob inequalities we have: 



E 



sup S r (x) 
re[o,t] 



< E 



sup (S r (x)Y 

re[0,i] 



< 2 I E 



(E [(5 t (x)) 2 ]) 



Therefore, we obtain the following bound for the expectation of A 3 (t): 



E 


sup A 3 (s) 






s£[0,t] 


V™ Jo 



(4.7) 



E 



(^(dv)^(dy) 7 (v + y) x t {0 . x] (v)) ds 



dx. 



Following the value of x we use different bounds: 
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On the one hand, for x < 1 we have t( ^(v) < (j) 2A £ and using the bound (v + y) x v 2X e < 
v 2X ~ £ y x , we obtain: 



„A-1 



3 



(tf(dv)tf{dy), (v + y) x t (Q . x] (v))d S 



dx 



< I x^ 1 {E 



^(dv)^(dy),^^) 



dx 



; 2 dx 



(4.8) 



o 



M A (/i?)M 2A _ e (/i?)ds 



On the other hand, for x > 1 we have l(o.a;]( w ) < ( f ) and using the bound (v + y) x v x < v x y x , 
we obtain: 



x^ 1 i E 



(^(dv)^(dy) 7 (v + y) x t^ x] (v))d S 



dx 



< 



x^ 1 i E 



+ 00 



x^~ l dx <E 



(4.9) 



i 
2 

|A| 



tf(dv)tf(dy), ^) 



{^ s {dv)^{dy),v x y x )ds 



dx 



Then, writing the right-hand side integral of (4.7) as the sum of the integrals on x <E (0, 1] and 
x e (1, +oo), gathering (4.8) and (4.9), we get 



E 



sup A 3 (s) 
«e[o,t] 



< 



/«o I 1 



M x {ii n s )M 2X _ e {ii n s )di 



(4.10) 



|A| 



[Mx(^)f ds 



Conclusion. 
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Gathering (4.2), (4.5), (4.6) and (4.10), we have: 



E 



sup d x (fi" 7 fi s ) 
se[o,t] 



< E 



d\(fj,Q , fio) + sup Ai(s)+ sup A 2 (s) + sup A 3 (s) 
se[o,t] se[o,t] se[o,t] 



2 A « /"* 
n|A| J 

8y^ f 1 

\/n I e 



E[M 2A ( M ™)] 



+ 7 



M A ( M ?)M 2A _ e (^)d a 



[M A (^)] 2 



According to Proposition A. 4 -(a), M Q! (/x" + /it) < M a (fiQ + fio) a.s. for any a € (— oo,0). Since 
Ho is deterministic, we get: 



E 



sup d\(fj% , fit) 
»e[o,t] 



+ TTT ) ) / 1 



(4.11) 



n|A| V n 



i(M A (/tf)M 2A _ e 0tf))* + ijM A (^) 



f 3. 



Finally, since Va6 < a + b and since M 2A (Mo) — A^a(Mo) + -^2A-s(Mo); we use tne Gronwall lemma 
to obtain 



E 



sup dx(f4,[i t ) 
te[o,T] 



< 



d x (n%, mo) + ^= M A (/tf ) + % M 2X _M) 



(4.12) x cxp 

where C x = + ^\^VT^ and C 2 = ^f 1 + f y/Tk 

This concludes the proof of Theorem 3.1 under (2.2). 



Kl 



T[ko + j^)M x (po 1 +Vo) 



5. Positive Case 

In the whole section, we assume that K satisfies (2.3) for some fixed A € (0, 1]. We fix e > 0, and 
we assume that Mo & -Mo H M.t l+e where 7 = max{2A, 4A — 1}. We denote by (fit)t>o the unique 
(fio, K, A)-weak solution to the Smoluchowski equation. We also consider the (n, K, Mo )-Marcus 
Lushnikov process, for some given initial condition fig = — $3 i=1 S Xi . 

We assume without loss of generality, for A e (0, 1/2), that e < h — A. Indeed, if e > | — A, it 
suffices to consider e < | — A, to apply Theorem 3.1 with e, and to use the bound M 2A+1 ?(/io + Mo) < 
M (fio + Mo) + M 2 \+ £ (fio + Mo) to conclude. 

We first present a lemma of which the proof is developed in the appendix. 
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Lemma 5.1. 

(5.1) 

Then, 



We introduce, for x £ (0, +oo) ; the following function: 



>(x) 



-7=l(0,i]O) + 



-2A- 



l(i,+oo)(a;)- 



(a) f +oo x^e» x) d x < 



(Hi) for (v,y) G (0, +oo) 

,,A 



A+e 
Xe y/n' 



+oo X X-1 



+y) ^ < 



- ^2A+ £ + 1 ^ [(w A y) 2X {v v ^A+e-^^ + („ A y){v V y)^-l lAg[1A1] ] . 

We set E n (t, x) = (x)~F^ (x) as defined in (2.6), for x G (0, +oo). We take the test function 
cj)(v) = l( x , +00 )(v). Since 

su Pd>o (i+Jj>> — (1 + %) A < +oo, we deduce that <fi G Ha- Again, 
computing the difference between equations (2.11) and (2.8) and using a symmetry argument for 
the first integral, we get 



E n (t,x) 



(5.2) 



E n (0,x) 



1 



2rc 



Ji<] 



(« - Ms) (*>) (m" + Ms) (dtf), (i4l (x , +00 )) («, </)}<te 

(Ai"(dw), (^l^.+oo)) (w,w)-ftT(w,w))ds 
/•+°° i / \ 



{j<JV(s-)> 

J(ds, <i(i, j), <iz). 



According to Lemma A. 2, we can write the first integral as: 

ft />+oo />+oo 

K(v, y) (Al (x>+oo) ) (v, y) (/i™ - Ms) (dv) + fi s ) (dy)ds 



JO 



n+oo />+oo 







t x>y K(x - y,y)t( Xi+00 )(v + y)- K(x, y)l(x, +00 )(v) 
+ I d x K(z,y)(At ) (z, V) dz\ (j% - fx s ) (dv) (/£ + n s ) (dy)ds 



o 



K(x - y,y) 

t p-\-oo 



l x>y 



1 



(x— y,+oo 



)(«) «-A*.)(d«) 



(p r s l + H S ) (dy)ds 



1(!E,+00)(«) (Ms - Ms) (d«) 



/o ./o 

pi />+oo />+oo 

d x K(z,y) (At {Xt+oo} ) (z,y) 



o Jo 



!(*,+«,) (u) (m™ ~ M«) 



(fx™ + fi s ) (dy)ds 



dz (fj% + fjt s ) (dy)ds. 
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Recalling that E n (s, x) = J + °° ^-( x ,+<x>)( v ) (m™ ~ Ms) (dv), we deduce that, 

I ft _ _ 

E n (t,x) = E n (0,x) + -J \B 1 (s,x) + B 2 (s,x)+B 3 (s,x)]ds 

(5-3) + J q jf^_ jf \ {At {x . +oo) ) (xi_,Xl_) I | „ jA x ; | { a - I 



where: 



J(ds,d(i,j),dz), 



r+oo 

Bx(s,x) = / [l^lf (a; - y,y)E n (s,x- y) - £„(s,a;)2f O™ + /z 8 ) (*/), 
Jo 

S 2 (s,x) = / / d x K(z,y) (At {x , +oc) ) (z,y)E n (s,z)dz + fi s ) {dy), 
Jo Jo 

— 1 f + °° 

B 3 (s,x) = — / K{v,v)[l {Xi+oo) (2v)-21 (Xi+oo) (v)]^(dv). 
n Jo 

Now, we apply the Ito formula to ipg(E n (t, x)), where ipe(-) £ C 2 (R) is an approximation of the 
absolute value function | • |. This function is chosen in such a way that: 

ip e (u) = \u\ if \u\ > 9; \u\ < tp g (u) < \u\ + 6 Vu E M: 

\ip' e {u)\ < 1 Vwel: sgn(uip' e {u)) = lVue R„; 

We(u)\ < \t { \u\<8} V«ei. 



(5.4) 



Furthermore, we consider for 6 the function defined by (5.1). We fix x G (0, +oo) and apply the 
Ito formula to ipg™ (E n (t,x)) (see for exemple [9]), 



9 nJE n (t, X )) = Lp e nJE n (0,x)) 

1 '* 



(5.5) +- / \B 1 (s,x) + B 2 {s,x) + B 3 (s,x)](p' d n (E n (s,x))ds 

* Jo (x) 



+M(t,x) + B±(t,x), 



where 



(p'e" x) (En(s-,x))j{ds,d(i,j),dz), 

-\ (At (Xi+oo) ) (.v; ..v: ) ^ (E n (s- x))|i j i { j<iv( s -)} 

J(ds, <i(«, j), <iz). 

Observe that, for all a; > 0, M(t,x) is a martingale whose expectation is equal to zero. 



16 



EDUARDO CEPEDA AND NICOLAS FOURNIER 



Now. we study the ^"-approximation of d\(fi , t l , fi t ): J Q + °° x x 1 ipg^ (E n (t,x))dx. According to 



(5.4) and Lemma 5.1 -(i), we have: 



(5.6) 



d\Qi",» s ) < / a: A ~Vfl» (E n (s,x))dx < d x (^,» s ) + 



Consider (5.5), integrate each term against x dx on (0 + oo), take the expectation: 

o+OO 



x A_1 IE 



(pg? (E n (t,x)) 



dx 



(5.7) 
where 



x x - V 9 » , (E n (0, x))dx + E [Bt(t) + B 2 (t) + B 3 (t) + B 4 (t)} 



Bi(t) = 

Bitf) = 

Bs(t) = 

B A (t) = 



1 
2 
1 
2 
1 

2 Jo 



+ OG ft 



^0 

oo />i 



x x 1 Bi(s,x)f'g n (E n (s, x)) ds dx, 



x A 1 B2(s,x)ip'g n (E n (s,x)) ds dx, 



o Jo 

OO pt 



x x 1 B^{s,x)(p'g n (E n (s, x)) ds dx, 



-foo 



x x ~ 1 B 4 (t,x)dx. 



We now study each term separately. 
Term B^t). 

We use the Fubini theorem to obtain: 

«-foo 



Bx{t) 



n+oo 



^x> y x x 1 (p' e „ {E n {s,x)) E n (s,x - y)K(x - y,y)dx 



H-oo 



x x 1 ip'gy. (E n (s,x)) E n (s,x)K(x,y) dx 



(i 



(fj% + ti s ) (dy)ds. 



Recalling (5.4), we immediately deduce that ip' gn (E n (s,x)) E n (s,x — y) < \E n (s, x — y)\, and 

(3=) 

ip'g„ (E n (s,x)) E n (s,x) = (p' e „ (E n {s,x)) \E n (s,x)\. Therefore, using the change of variable 
x t— > u + y in the first integral, we get: 



Bi(t) < 



1 



r 


" /-+oc 




.Jo 


/• + 00 






A-l 




•70 




/•-t-oo 


/•-t-oo 



(n + y)*- 1 |K(s, u) I K"(u, y)du 
ip' en (E n (s,x)) \E n (s,x)\K(x,y)dx 



Jo 



|£„(a,z)| (z + y)^ 1 - ^' (£„(s,z)) 



.A-l 
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Recall again (5.4). Since \E n (s,z)\ > 9?-. implies ip' e „ (E n (s,z)) = 1, and since (z + y) 
z x ~ x < 0, 



A-l 



\E n (s,z)\ [{z + y) x - 1 -\^ %) {E n (s,z))\z x - l \ < \E n (s, z)\(z + y) x ~H z)|<e „ J 



< 



A-l 



Therefore, using (2.3): 



*(*) < f 



< H 

~ 2 



'0 Jo 

ft /.+oo 



qn ,2A-1, 



r [z) z x -Uz 



(fj,s+H s ) (dy)ds. 



) z ZA - l dz + y x 

'o Jo 

We used (z + y) 2A_1 = (z + y) A (z + 2/) A_1 < (-z A + y A ) z . Finally, according to Lemma h.l-(i) 
and (mJ, we get: 

ft />+oo [( 



Bi(t) < 



(5.8) 

Term B 2 (t). 



< 



K0 

2 Jo Jo 
kq(X + e) 



^±|1(i + 2/ a ) ( M » + Ma )(iy)d» 

[Mo(m? + M-) + Ma(m? + M-)] d«- 



First, observe that 

|(^l(x,+oo)) (Z,V)\ = |l(x,+oo)(^ + V) - l(x.+co)(z) - l(x,+»)(y)| 

(5-9) = l{xe(0,zAy)} + l{a-e(zVi/,2 + y)} j 

whence, 



x x 1 1(^4.1(^+00)) (z,y)\dx 



zAy 



<- x ~ l dx 



z+y 



zVy 



~ x - l dx 



(5.10) 



< j{zhy) x . 



Thus, recalling (5.4), we get: 

rt />+oo />+oo 



5a(*) < = 



< 



Ki 

A 

Kl 

A 



Jo Jo 

• t />+oo />+oo 



y)| [ - (2 A y) ) (tf + fi s )(dy)dzds 



\E n (s, z)\z x ~ 1 y x {pi n s + n„) (dy)dz ds 
c2a(aC Us) M x {^ + ^ s )ds. 



JO Jo 

t 



(5.11) < 

We used (2.3). 
Term B 3 (t). 

Remark that | (Al( Xt+00 )) (v,v)\ = \l( Xi+oo) (2v) - 21( Xi+oo) (v)\ < !{„>§}• 
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Since J +O ° t{ v> *}X X 1 dx 



A 



-, we deduce: 

f t />+oo 



Bs(*) < ^ J x ^ J J K(v,v)\(Al {X!+oo) )(v,v)\fj%(dv)dsdx 



(5.12) 



< 



< 



2Xn 



ds 
o 



(2») 2 V?(d«) 
M 2X {^)ds. 



We used (5.4) and (2.3). 
Term £ 4 (i). 

First, remark that from (5.4) we have ip'gn (z) 
Lagrange inequality, 

<pgn m) (E n (s,x) + i (Al (a , +oo) ) (.VI..\ ; i)j - ^jE n (s,x)) 



< for all z, whence, due to the Taylor- 

0») 



< 



'(x) 



(Ai (:c , +00) )(x:,xf) 



Then, 



E[B 4 (f)] < 



+ 00 



x x ^E 



+oo 2 
an 

10 Ji<j JO a ( x ) 
1 



- (^(x.+Oo)) (^s_,^_) 



1 



{j<N(s-)} 



z< 



«r(x»_,xf_) I 



J(ds, d(i,j),dz) 



d.i: 



< 



< 



2 

n 

2ko 



E 



o 



^ E 

i<j<N(s) 



K{ n2 l gn Xi) [(Al (x , +oo) )(XlXi)] 2 dx 

(as) 



E 



- E (*:+^) ; 



( K ) i<j<N(s) 



{^x<Xi^,X 3 s + ■"■X*VJf|<x<X*+Xi 



(/;)' 



We used (2.3) and (5.9) (since the sets are disjoint, the product of indicators vanishes). Therefore, 
using that (v + y) x < v x + y x and a symmetry argument, we get 



ELB 4 (*)] < — / E 
n 



rf(dv)tf(dy)y 



'(x) 



(1 



v-\-y) dx 



ds. 
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According to Lemma 5.1-(iii), and since (v A y) a (v V y)? < v a y@ + y a v 13 for a > and f3 > 0, we 
have 



'(a:) 



(lx<«Aj/ + l»Vy<a;<w+y) dx ) < 



<M?(d«K(dy),«V) 



Finally, we deduce the bound: 



E[B 4 (t)] < 



8k 



[M x (tf)f + C [MM)M 2X+£ W)] l Ag (o,i/2) 



(5.13) 



-C [Mi(/x")M 4 A+e_l(/X™)] 1ag[1/2,1] 



ds, 



where C* = (A2 2A+e + l). 
Conclusion. 

Gathering (5.8), (5.11), (5.12) and (5.13), from (5.7), we get: 

«+oo 



E[dA(/i?,/i*)] < 
, K (X + e) 



Xe^/ri 

^ ftp 

nA 
8C^q /* 



x A -Ve" ) (-Er l (0,x))da; 



■[A/ (aC + M.) + A/ a (m" + Vs)]ds ■ 



A 7 



E[d A «,// fl )Af A (^+// s )]ds 



E[M 2A (/0] ds 



8k 
X^/n 



E[M x (^)fds 



Xy/n Jo 



[MM) M 2X+e (tf)] l A6(0ll /2) + [MxCOM^-iCO] 1ag[i/2,i] 



(is. 



We use (5.6) to bound the first term on the right-hand side. According to Proposition A. 4 -(a), 
M a (^ + fJ. s ) < M a (fiQ +/io) a.s. for a < 1. Since /Xq is deterministic, we get (recall that 2A + e < 1 
if A e (0,1/2)): 

E [d A (M?, A*t)] < <M, Mo) + t—i= + tK °} X t_ £) (M (^ + Mo ) + M x QjQ + Mo)) 



«i 



Ay 7 ?! Xey/n 
t 2 2A_1 



+yM A K + Mo ) / E [d A (/C Ms)] ^ + 



'HI 



hA 



E[M 2A (^)]rf s + -^[MA(^)] 2 
Av« 



AVh 

+ ^|Mi(^) / E[M 4 A +£ -l(M?)]l AG [1/2,1]^. 



[M 2 A(Aio)M2A+ £ (Aio)]lAG(0,l/2) 
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Again, according to Proposition A. 4 -(b), E [M a (/x™)] < M q (mo ) exp[s (7x,aMx(Mo)] for a > 1, and 
where Cx, a is a positive constant depending on A, a and k . Thus 



Recall that 7 = max{2A,4A - 1}. Observe that for fi £ M + , M a (ji) < M (fi) + Mp(p) for any 
< a < p. Elementary computations allow us to get: 



for some positive constant Cx,e depending on A, e, kq and K\. We conclude using the Gronwall 
lemma that Theorem 3.1 holds under (2.3). 



6. Special Case 

Now we are going to study the special case (2.4) for which A £ (0, 1]. We have a better result 
and a simpler proof than (2.3). 

In the whole section, we assume that K satisfies (2.4) for some fixed A £ (0, 1]. We fix e > 0, and 
we assume that mo £ A4^ fl M. J\+ e - We denote by (Mt)*>o the unique (mo > K, A)-weak solution to 
the Smoluchowski equation. We also consider the (n, K, Mo )-Marcus Lushnikov process, for some 
given initial condition Mo = i ^* ■ 

As we did before we introduce E n {t, x) = F Mt (x) — F^ 1 (x) for x £ (0, +00), as defined in (2.6). 
We observe that l^+oo) <= H\, since sup tI>0 v~ X \^( x ,+oo){ v )\ = x < +00. Exactly as in Section 
5 (see (5.3), take the absolute value and integrate against x x ~~ 1 dx), we obtain: 



E [dx(tf, ih)] < dx(^, mo) + + tK °} X t- £) (M„(/tf + no) + M x (f% + Mo)) 




Jo 




E [d x (jj?,iH)] < dx(^,Vo) + (1 + 1)-^ (l + [Mo(Mo + Mo)] 2 + [M 7+e (Mo + Mo)] 2 ) 
x cxp [t C x , s M x {Ho + Mo)] + Ca, e M a (mo + Mo) / E [d x (^, ft B )] da, 



(6.1) 



dx ((4, Mi ) < dx (Mo , Mo) + Cx (t) + C 2 (t) + C 3 (t) + C 4 (t) , 
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where 
Ci(t) 

c 3 (t) 
c 4 (t) 



t />+oo />+oo 



Jo Jo 



„A-1 



l a: > 1/ iir(x - y,y)\E n (s,x - y)\ + \E n (s,x)\K{x,y) 



dx 



t f-hoo /*+oo /"+00 



JO JO Jo 



« + fj, s ) (dy) ds, 
x x ~ 1 \d x K{z,y)\ \ (At (Xt+oo) ) (z,y)\ \E n (s, z)\dz dx 

(/£ + M.) (<*ff) da, 



i /*+CO /*+G© 



2rc 



o Jo Jo 



x x 1 K(v,v) |l(a. )+00 )(2u) - 21( X)+00 )(«)| dx n™(dv)ds, 



./i<j ^0 



J(ds, d(i,j), dz) 



dx. 



We now study each term separately. 
Term d{t). 

We have, using the change of variable x H> w + y, (2.4) and using the fact that x A_1 is 
non-increasing function: 



Ci(t) < 



< 



/- + 00 /■+(» 

4 / (^ + y) A - 1 (^Ay) A |£; n ( S ,u;)|d W (/i:+^)(d 2 /)d S 

z Jo Jo Jo 

ft p+co r+oo 

+ 4 / ^- 1 ^/\y) X \En(s,x)\dx(^ + fi s )(dy)ds 

1 Jo Jo Jo 

/ «/- VK(s, w)\ dw (/x" + Ms ) (dy) ds 
Jo 

(.+oo /■+(» 

+ ^ / / / x x - 1 y x \E n (s,x)\dx(^ +^ s )(dy)ds 
1 Jo Jo Jo 



< k q / M x (n™ + fj, s )d x (n™,Vs)ds. 
Jo 



(6.2) 

Term C 2 (t). 

Recall (5.10), use (2.4), we have immediately: 



t p+oo />+oo 



c 2 (t) < - 



\d x K(z, y)\ (z A y) X \E n (s, z)\ (fi + Ma ) (dy) 



'0 Jo JO 

• i z'+oo /'-(-oo 



< 



«1 

Jo Jo 



|£7 n (s, z)\z x ~ l y x + n.) (dy)dz ds 



(6.3) 

Term C 3 {t): 



d x (fJ%, fJ, s ) Mx((J% + (J. s )ds. 
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As before, recalling (5.12), we write: 
(6.4) C 3 (t) < '- 

Term C 4 (i): 



A n 



M 2X (^)ds. 



The submartingale term is going to be treated exactly as in the case A < 0. Using similar 
arguments as for the term A^t), we get 



E 



sup Ca(s) 

s6[0,t] 



< 



„A-1 



E 



V%{dv)ii™{dy) , K(y,y) [(At 



dx. 



Using now (5.9) and (2.4), we deduce that 

4^ 



E 



sup 6*4(5) 

sG[0,t] 



< 

















E 


[I 









^(dv)^(dy),(v/\yy 



(6.5) 



First assume that x < 1. Since ^{ x e(o,vAy)} ^ 



[l{a;e(0,t)Ay)} + ^{x£(Wy,v+y)}] ) 

(vAy) x 



dx. 



, since t^ x£ ^ vVyv+y ^ — 



(v + y) } 



< 



(v V y) x 

2- t , and since (v A y) x (v A y) x < v x y x and (w A y) x (v V y) A = v x y x 7 we deduce that 



(^(dv)fi"{dy) , (u A y) A [l {xe(0 ^ Aj/ )} + t {xe{vVy ^ v+y)} ] ) 



Thus, 



„A-1 



E 



{^{dv)p n s {dy) , (vAy) A [l {a . e(0it , Ay)} + l{ a:e (wj / ,«+y)}] ) 



dx 



(6.6) 



< Vl + 2 A / x^~ x dx 



[Mx{n n s )] 2 ds 



2Vl + 2 A 
A 



[M A (^)]^s 



(wAy) 2A+£ (w + y) 2A+£ 

Next consider a: > 1. Since 1{ x& (o,va v )} < — 2 \+e — , and l{ xe ( vVy ,v+y)} < — ^ 2A+e — < 

(v V y) 2A+£ 

2 2A+e V ^ and since ^ A y j\( v A y )2A+e < w A y 2A+e and ^ A y )A( w y y ^2\+e < y \ y 2\+s + 



r 2A+£ 



v 2X+£ y x , and using the symmetry, we deduce that 
{n™{dv)n n s {dy) , (v A y) A [l {:re(0it , Aa)} + t {x<£[vVytV+y)} ] ) 



< 



;i + 2 2A+£ + 1 ) 

T 2X+e 



M x (tf)M 2X+ M. 
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Thus, 



/>+oo 



„A-1 



E 



(fi%(dv)(J%(dy) , (wAy) A [1^(0,^)}+ t {xe(vV y, v+y)} ]) dt 



d;r 



(6.7) 



< \/l + 2 2A + e +! 
2>/l + 2 2A + £ +! 



f+oo 



x 2 da; x I E 



M A (/^)M 2A+e (/^)ds 



E 



M A (/£)M!u +e 0£)da 



Gathering (6.5), (6.6) and (6.7), we obtain: 

'ko f Vl + 2 A 



S 



SUp C^s) 
s6[0,t] 



< 



A 



E 



(6.8) 



^ 1 + 2 2A+e+l 



s 



M x (tf)M 2X+ Mds 



Conclusion. 

Therefore, gathering (6.2), (6.3), (6.4) and (6.8), we obtain: 



E 



sup dx(n^,n s ) 

sG[0,t] 



< 



„t 



fn 1 A 



E[M 2A (/*")] 



E 



+ VTT^tt (e 



[M A «)] 2 ds 



Af A (// s l )M 2A+£ (// s l )d S 



Observe that M a (// l ) < M (/x l ) + M 2A+e ( / U l ) for a = A, 2A. Proposition A.4 implies that 
for a € (0,1], M a (fi? + Mt) < M^ft + Mo) a.s. and for a = 2A, 2A + e, E[Af Q (/x?)] < 
A/q^/Xq) exp[s C A . Q Af A (/x 1 )] where C A;Q is a positive constant depending on A, a, «o and K\. Since 
/io is deterministic, we deduce that 



E 



sup d A (/i" , Us) 
»e[o,t] 



< 



d A (^, mo) + (1 + 1)% (M (/#) + M 2A+e ( M ^)) exp[*C A , e M A 

-v / 77, 



for some positive constant C AjE depending on A, e, ko an d Ki- We conclude using the Gronwall 
lemma. 



7. Choice of the initial condition 

The aim of this section is to prove Proposition 3.2. We thus fix A e (— oo, 1] \ {0} and [Iq E 
.M A n M^x- We first treat the case where /xq is atomless, next the case where /xq is discrete. 
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7.1. Continuum System. — We assume that /io is atomless. For < a < A < +oo, we consider 
/j,o\k , the restriction of /io to K = [a, A}. We consider also TV points a = x < xi < ■ ■ ■ < xjy < A 
such that: 



(7.1) fj, ([xi-i, Xi)) = —, V i = 1, ■■■ ,N and fj, ([x N , A]) < 



1 



We will use the points {xi}i = i.... . n to construct the discrete measure /ig choosing a and A following 
the value of A function of n. 

7.1.1. Case A £ (— oo,0): — First, we choose a n < A n as follows: 

1 



(7.2) 



and 



x A no(dx) < —=. 



Next, we assign the weight /io %i)) = h to the point Xi and we set 



= ^ 



i=l 



(7.3) 

If a < 0, we get: 

Mating) = -y^xf = Vxf Hodxi-x, Xi)) = V / if l [a!i _ lia! .)(a;)^o(da;) 
n i= i i=i 

(7-4) < / * a l[xi-i,x 4 )(aOA*o(<&0 



i=l 



x a /j, (dx) < M a (fj, Q ). 



For the distance, we hcLve, with ^ n — [*i n ,-A n ]; 



dAGuok^Mo) 



«A-1 



da; 



+oo 



X [/i ((0, x)) l x <a„ + /^O ((^n, X)) lz>A„ + jL*0 ((0, 0„)) l x >aj dx 
H-oo />+oo 

x A_1 dx Ho{dy) 



o J 2/ 

a n /*+ao 



< 2 



(7.5) 



< 



2a^ 
|A| 



x A 1 dxn (dy) 

x A_1 dx Ho(dy) 

x^ 1 dx(j> (dy) + 
f + °° 1 /" 



Ja„ 
a« /■+00 







+oo />+oo 

/ x x ~ x dx /j, (dy) 

A n Jy 
+oo 1 

X, ' • - - 1 



y A Mo (dy) < jYj 7= (2M 2A (/i ) + 1) , 



we used (7.2) for the last inequality. Next, we introduce the notation i x = max{i : Xi < x; i = 
0, • • • , N n } for x > a n . We remark that Hq ((0, x]) = if x < a„ and /Zq ((0, x]) = (j,q ((a n , XjJ) if 
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x > a n . Hence, 



^A(Mo>Mok„) = / x 



< 



+00 



„A-1 



+00 



^(o,x)(y) (Mo - MokJ (dj/) 



A„ 



(7.6) 



a; I /■»<)( fan, a*,)) - fJ-o([a n ,x))\ dx 
x A_1 |Mo(fan,£iJ) - ^o([a n , Ai))\ dx 



A n r+oo 

x x ~ 1 ii Q {{x ix ,x))dx + / x x ~ 1 ii {[x Nn ,A n ))dx 

2 f + °° \ i 2 \ 2 

nj„ An " AK/n 



We used |/Zo(fan, 3%,)) - /i ([a„,a;))| = ^((i;,,!)) < nod^-i,^')) < ^ for some j = 1, 
and (7.2). Finally, from (7.5) and (7.6), we obtain: 



dx(fJ,Q , fio) < ^(^ojMoIO + ^(MokniMo) < mV (2M 2A (mo) + 3) 



7J.£. Case A € (0, 1]: — First, we choose a n < A n as follows: 



(7.7) 



x x /io(dx) < —= 



and 



A n = (y/n) ' 



Next, we assign the weight /Xo([^i-i, Xij) = — to the point Xi-i, recall that xq = a n . We set 



W„-l 



(7.8) 



If a > 0, we get: 



H™(dx) = - V s x 

11 — J 



i=0 



N n -1 N„ 



N„ 



M Q (j#) = - ]T x? = ^xf^ fi ([xi-i, Xi)) = ]T I x°'_ 1 t [Xi _ uXi) {x)^ (dx) 

i=0 t=l i=l 

"n p+co f x N n 

(7.9) < / ^l^^oWwW = / x>o(^) < M a (/i ). 
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For the distance, we have, with K n = [a„, A n \: 



d\(fJ,o\K n , Mo) 



„A-1 



l[:c,+oo)(y) (Mo|a„ - Mo) (dy) 



dx 



+00 



X [tt ([x, O n )) lx<a n + MO ([a?, +00)) 1 K >A„ + MO ([Ai, +°°)) 1U<aJ 



x A 1 no(dy) + 



JO 

+00 pA n 



+00 />y 



A n JA„ 



x A 1 dx fia(dy) 



x dx fio(dy) 



< 



A„ JO 
a n ry 



JO 



x A " 1 dxn a (dy) + 2 



+00 r-y 



x x 1 dx no(dy) 







(7.10) 



1 

< - 
" A 



x X jio(dx) 



2A~ X f +oa 1 

y 2X Mdy) < V7 =(l + 2M 2A (/i )) ) 



A 







we used (7.7) for the last inequality. Next, using the notation i x = min{i : Xi > x\ i = 0, • • • , N— 1} 
for x > a n , we remark that ([x, +oo)) = if x > A n and /Xp ([#, +oo)) = /xo ([^i x , A n )) if a; < A n . 
Hence, 



^(Mo^oIa,, 



-foe 



i(x,+oo)(y) (mo - MokJ (dy) 



(Z.r 



x \no((x ia: ,A n )) - fi ((x,A n ))\dx 



+ x x 1 \fio([xi x , A n )) - /x ((a„,A„))| 
Jo 



(7.11) 



x ii ((x,x ix ))dx < 



x x ~ 1 dx=-^A x <-± r , 
An A 



we used \^o{{x ia: , A n )) - fj, ((x, A n ))\ = /i ((x,x ?; J) < (i ([xj-i, x 3 )) for some j = 1, • • • , N, and 
(7.7). Finally, from (7.10) and (7.11), we deduce: 



d\(fJ-o,Vo) < 



(M 2A (/i ) + 1) . 



7.2. Discrete System. 

(7.12) 

We set for i£N: 
(7.13) 



Let us thus, consider zxq 6 A4 + with support in N, i.e. 



u Q = ^2 a k 6 k , with a k £ 

k>l 



Mo 



k=l 



7.2.1. Case A <G (— oo,0): — We choose A n such that: 



(7.14) 



fc>A„ 
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and we set, 
(7.15) 

(7.16) 



j An 

^ = - J2a n k S kl with 



fc=l 



a™ = L na iJ i 



\n(o>i H h a k )\ - [n(ai H h afc_i)J for fc = 2, • • • , A„, 



where [-J is the floor function. Remark that chosen in this way, the a k are non-negative integers 
and (j,q can be written as — y^,-!*] S X4 , hence is the measure we search. Observe that for 
& = !,••• , A„ , we have 



(7.17) 

If a < 0, we have: 



(a? + ..- + ajf)-(ai + "- + a fc ) 
(ai H h «fc)J - (aH ha*) 



< 



M a (n r 



1 A„ 
k=l 

1 1 A " 



fc=2 



— InaiJ H — V^lnfai 

r» n ^- — ^ 



fe=2 

c )Jfc Q -- 51 L"("i + --- + a fc )J(fc + i) c 

k=l 

= - ( [na x \ + AI [n(ai + • • • + a A JJ - 2 Q hn«ij ) 



fc=2 



1 



A„-l 



+- V Ln(ai + • • • + a fe )J (fc Q - (fc + 1)«) 

n * — » 



fc=2 



A„-l 



< oi (1 - 2 Q ) + A£ (ai + • • ■ + a A J + X! (ai + • • ■ + a k ) {k a - (fc + 1) Q ) 



*;=2 



A„-l 

fe=i 



(7.18) 

Next, for the distance, we have: 



A n -1 

4" + £ CT - C? + !)") 



A„ 



fc=l 



< 



1(0,*) (y) (/"o 4 " - Mo) (dy) dx 

/ xA_1 / E a kS k {dy)dx = V a fe / o^cfec 
Jo Jo fe>An 4 Jfc 



(7.19) 



i^E^ A ^m^> 



|AlVn' 
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we used (7.14). Next, 



+00 



An-1 r k+l 



= E 



fe=i 



„A-1 



l(o,)(lf)(l<S-/'o-)(*) 



dx 



i=l 



+00 

dx + I x*' 1 

A„ 



dx 



2 r+00 2 

< - / x^dx < —-, 
nj 1 \X\n 



(7.20) 

we used (7.17) for the last inequality. Finally, from (7.19) and (7.20), we have: 

dx (Mo>/«o) < d x (mo,^") + (^ n ,Mo) < fl + ^= 

7.2.2. Case A G (0, 1]: — We set 



1. 



(7.21) A n 
Note that chosen in this way, we have A~ x < implying 

(7.22) V a k k x < A- x V a k k 2X < -Lm 2A ( Mo ). 

k>A n k>A n v 

We set the measure Hq as defined in (7.15), with 

for k = 1 , • • • , A n . 



(7.23) 


«2 = 


n a. 




n ^ eg 






i>k 




i>k+l 



Observe that, since J2k>i a k = Mo(p-o) < ^h(fJ-o) = J2k>i a kk X < +00, the weights {a k l } k >i are 
well-defined. Remark that chosen in this way, the a k are non-negative integers and can be 
written as ^ X^Ji ^i, hence /Uq is the measure we search. 



For 1 < j < A n , we have: 



fc=j 



-a k - a k 





1 

n 


n E ai 




i>j 




1 

n 


A n 


< 






i=j 







A„ 


- E 


a. 


~E afc 


i>A„+ 


1 


k=j 


A„ 






-E afc 













(7.24) 



< - 



A„ 



A„ 

E afc 



2 

< -. 
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If a > 0, we have: 
MM) = 



, A„ A n 
n ^ K n ^ 



k=l 
-1, 



fe=l 



i>k 



1 A„ 

n ^ 

k=l 



i>k+l 



k=2 



i>k 



[k a — (k — l) a ] + 



n ^2 a i 

i>i 



4" 
n 



k a 



n eg 

i>A n + l 



(7.25) 



fc>l \i>k 

For the distance, we have : 



i>i k=i 



(mqSmo) 



[a;,+oo) 



(y) (mo " - Mo) (dy) 

/ xA " 1 / X! a k 6 k (dy)dx = a k / * A_1 



fc>A„ 



(7.26) 

we used (7.22). Next, 



k>A T 







r+oo 


/•+00 


f n A„\ 












/ ^ 








Jo 


Jo 



i [*,+«>) (y) (mo - Mo 4 ") ( rf y) 



c/.r 



= X 

3=1 



.A-l 



(7.27) 



< 



j'-i 



-a k - a k 



dx 



x x ~ x dx < 



2A* 



< 



o 



A n A \/n 



we used (7.24) and (7.21). Finally, from (7.26) and (7.27), we obtain: 

^a(Mo> Mo) < d^,4») + d x (4'Wo) < ^=(M 2A ( Mo )+4). 

7.3. Conclusion. — In any case, (A e (-co, 1] \ {0} and fio either atomless or with support in 
N), we have built a measure of the form //q = i Xa=Ti satisfying the desired conditions on the 
moments and distance. It is straightforward to show that N n = n (^(dx), 1). Hence, according 
to (7.4), (7.9), (7.18) and (7.25), we deduce, 

(7.28) N n = nMo(tfS) < nM (p ). 

This concludes the proof of Proposition 3.2. 



Appendix A 

This section is devoted to some technical issues. 
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Lemma A.l. — Consider A G (— oo, 1] \ {0}. Then, there exists a positive constant C^, depending 
on <j) and A such that 

( i/Ae(-oo,0), (x + y) x \(A^)(x,y)\ < C^(xy) x V^%, 
(A.l) { i/Ae(0,l], (.T + y) A |(^)(x,y)| <C^(l + x 2X +y 2X ) V0 G H\, 

{ if \€ (0,1], (xAy) x \(A(j))(x,y)\ < C r j ) (xy) x Vct>eH e x . 

Proof. — Assume first that A G (— oo,0) and 4> G rl\. Since \4>{x)\ < Cx x for some constant 
C > 0, we have 

[x + y) x \{Act>){x,y)\ <C(x x Ay x ) [(x + y) x + x x + y x ] <C(xy) x . 
Next, for A G (0, 1] and <p € rl\, since \<fi(x)\ < (7(1 + x x ) for some constant C > 0, we have 

(x + y) x \{A<j>){x, y)\ < C(x + y) x [3 + (x + y) x + x x + y x ] < C (l + x 2X + y 2X ) . 
Finally, for A G (0, 1] and <f>(x) G H\, there exists C > such that \4>(x)\ < Cx x and we have 
(x A 2/) A |(A0)(x, y)\ < C(x A y) A [(x + y) x + x x + y x ] < C(xy) x . 

□ 

Lemma A. 2. — Let A G (-oo, 1] \ {0} and K G VF 1 ' 00 ((e, 1/e) 2 ) for every e G (0,1). If K 
satisfies (2.2), then for all (x,v,y) G (0,+oo) 3 : 

K{v,y) [l( QtX ](v + y) - l(p lX ](v) - l(o,x](?/)] 
(A.2) = K(x-y,y)t i0tX] (v + y)-K(x,y)l( 0iX] (v) 

«+oo 

d x K(z,y) [l(o,x](2 + 2/) - l(o,x](z) - 1(0,1] (y)] dz. 

//A' satisfies (2.3) or (2.4), then for all (x,v,y) G (0,+oo) 3 : 

K(v,y) [l( x , +00 ){v + y) - t( x , +00 )(v) ~ t (x,+oo)(v)] 
(A.3) = K(x ~ y, y)lx>yl(x,+oo)(w +y)~ K(x, J/)l( x ,+oo)(«) 

+ / d x K(z,y) [l( x ,+o )(^ + y)-l(x,+oo)(^)-l (y)] <fe. 
Jo 

Proof. — For A G (— oo, 1] \ {0} we have that K(-, •) and its weak partial derivatives belong to 
L°° ((e, 1/e) 2 ), whence, for all < a < b < +oo and for all y > (see for exemple [16]): 

(A.4) f d x K(z, y)dz = K(b, y) - K(a, y). 



First assume (2.2), and fix A £ (— oo,0). Remark that: 

/•+00 pb 

/ d x K(z,y)dz= lim / d x K(z,y)dz = lim K(b,y) - K(a,y) = -K(a,y). 

J a b^+oo J a b^+oo 



Hence, 



+00 /'+00 

d x K(z, y)t( ,x] (z + y)dz = t x>y t v < x - y d x K(z,y)l v < z < x - y dz 

Jo 

= ^(o,x](v + y)[K(x-y,y)- K(v,y)]. 
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Next, 



+ 00 /"+00 

d x K(z,y)l(o !X ](z)dz = -t v < x / d x K(z,y)t v < z < x dz 

Jo 

= %(o,x](v)[K(y,y)-K(x,y)], 

+ OO /"+QO 

d x K(z,y)t^^(y)dz = -l y < x / d x K(z,y)l v < z dz 

Jo 

= %(o,x](y)K(v,y). 

Adding these three terms to the terms on the right-hand of (A. 2) the result follows. 
Next, assume (2.3) or (2.4). Observe that for (x,y, z) G (0,+oo) 3 , we have 

(A. 5) ~^z>x — y ~^-y>x — ~^-y<x^- Z >x— y 

Thus, 



/ d x K(z, y) [l z + y >x - t z>x - l y>x ] dz 
Jo 

= / d x K(z,y)[l y < x l z>x -y -t z>x ]dz 
Jo 



= l y < x l v>x -y / d x K(z,y)dz - t v>x / d x K(z,y)dz 

J x—y J x 

= l v < x t v>x -y [K(v, y) - K(x - y, y)] - t v>x [K(v, y) - K(x, y)] 
= [l v >x-y - ^ y >x] K(v, y) - t y<x t v>x - y K{x - y, y) 

-t v>x [K(v,y)-K(x,y)}. 
Adding these terms to the remaining terms on the right-hand of (A. 3), the result follows. □ 

Now we will show a lemma which is useful to show Proposition A. 4 stating that the a-moments of 
Ho and fi^ remain bounded in time. 

Lemma A. 3. — Consider a £ 1 J 6 (— 00, 1] and a kernel K satisfying either (2.2), (2.3) or 
(2.4). We set = x a . Then, 

(1) tfae(-oo,l], (A$)(x,y) < 0, for (x,y) e (0,+oo) 2 , 

(it) if a € (l,+oo), K(x,y)\(A0)(x,y)\ < C x , a (x a y x + x x y a ), for (x,y) e (0,+oo) 2 , 
where Cx, a is & positive constant depending on X, a and kq. 

Proof. — Point (i) is obvious, since for a < 1, (x + y) a — x a — y a < (x a + y a ) — x a — y a = 0. 
Next, if a > 1, using (2.2), (2.3) or (2.4), there holds K(x, y) < k q (x x + y x ). We get 

K(x,y)\(Aff)(x,y)\ < n x x [\{x + y) a ~ x a \ + y a ] + K y x [\(x + y) a - y a \ + x a ] 

< a Ka [(x x y a + x a y x ) + (x + y)"" 1 (x x y + xy x )] 

< C [{x x y a + x a y x ) + (x*- 1 + y^ 1 ) (x x y + xy x )] 

< C (x x y a +x a y x +x x+a - 1 y + xy x+a - 1 ) . 

Note that x x+a ~ 1 y = x a y x (| ) 1_A = x x y a f^j" < x a y x t x>y + x x y a t x < y . We have an equiv- 
alent bound for the fourth term and the result follows. □ 
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Proposition A. 4- — Consider A £ (— oo, 1] \ {0} and a coagulation kernel K satisfying either 
(2.2), (2.3) or (2.4). Let € and denote by {nt)te[o,T) the (/zo, K, \)-weak solution to 

Smoluchowski's equation. Let [1q be a deterministic discrete measure and (/i")t>o the associated 
(n, K, fiy)- Marcus- Lushnikov process. Leta<EM., then 

(a) if a < 1, in- M a (pt) o-nd t ^ M a (fi") are a.s. non-increasing; 

(b) if a > 1, there exists a positive constant C\_ a depending on X, a and kq such that M a (fi t ) < 
M a (jto)exp[tC x , a M x (jio)] andE[M a (n?)] < M a (/iff)exp[*C A , a M A (^)]. 

Proof. — Let (j>(x) = x a . For point (a), first consider (2.8). From Lemma A. 3.— (i), we immediately 
deduce 

^{lH(dx),<t>(x)) = = ^(vH(dx)iH(dy),(A<l>) (x,y)K{x,y)) < 0. 

Next, consider (2.10) and remark that <j> (x\_ + X J S _} - <p {X\_) - <j> (x j s \ = {A<t>) (xi_,X J 8 _^ . 

From Lemma A.3.-(i) and since J is a positive measure, we deduce that the jumps of M Q (/z") = 
(fj%(dx), 4>{x)) are negative and the conclusion follows. 

For point (b), consider (2.8). According to Lemma A. 3.-(ii), we deduce: 
j t M a ( Pt ) = \(nt(dx)iit(dy),(A0)(x,y)K(x,y)) < ^ (fx t (dx)fx t (dy), x a y x + x x y a ) 

< Cx, a M x (jit)M a (jM) 

< C\ >a M\ (/io) M a (pt) , 
we used the point (a). We conclude using the Gronwall lemma. 

Next, we take the expectation in (2.11). Remarking that (A<p) (x,x) > 0, using Lemma A. 3.-(ii), 
since //q is deterministic, and since M a (/it") = (fj,f(dx), 4>(x)), we deduce: 

1 '* 



E[M a (tf)] = M a (jj,o) + - I E[((j,™(dx)fi™(dy),(A<t>) (x,y)K(x,y))]ds 

1 r\ 



2n 



E [(${dx), (Acj)) (x, x)K(x, x))} ds, 



< M Q (^) + % [ E[(tf(dx)tf(dy),x<*y x + x x y>*)]ds 



< M a (jj$)+C x , a / K[M x (tf)M a (tf)]ds 



< M a (fi l ) + C x , a M\ (/#) / E [M a (/£)] ds, 

Jo 

where we used the point (a). We conclude using the Gronwall lemma. □ 
Finally, we present the 

Proof of Lemma 5.1. — Assume that A G (0,1] and recall (5.1). First, for (i) and (ii), by direct 
integration, we have 

x^Bl'dx = —= / x x ~ l dx + - F = / x- A_s_1 da! 

1 1 2 

< 



o 



Xy/n (A + e)y/ri \\fn~ 1 
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and 



o 



+00 1 Z" 1 1 r+ ao 

x 2A_1 0?^dx = -= / x 2X ~ 1 dx + -= / x'^dx 



1 , 1 (A + e) 



2Xy/n E\fn \e\fn 
Next, for (Hi) we have: 

r+°o r v ^v x x ~ x f v+y x A_1 

A n {v,y)=V X (t X<vA y + l V Vy <X<V + y)dX = V^J — ~ fife + U A / "^J" ^ 

JO ^(a) JO tf (a,) i»vy W (» 

:= I n {v,y) + J n (v,y). 

We have the following bounds: if v A y < 1, then 

/„(«, y) = V A f Vi^-^a; = ^V(« A 2/) A < ^Vy\ 
Jo A A 

Next, if v A y > 1, 

I n (v,y) = v x / y/nx^dx + v x / Vna 3A+e_1 dx 



< A/nu A 



1 (tiAjf) 



3A+e 



A 3X + E 



< ^ [v x + (vAy) 3X+s V X ] 



A 

< ^ [*V + (v A y) 2A ( V V y) 2A+£ l Ae (o,i/2) + (« A */)(« V y) 4A+£ - 1 l Ae[ i/ 2 ,i]] ■ 
Thus, in any case 

In(v, I/) = x ^ A + ( v A ^ 2A ( V v J/) 2A+£ lA e( o,i/2) + (« A y)(t> V y) 4A+£ - 1 l Ae [i/ 2 ,i]] ■ 
Next, since a; A_1 and 9? are non-increasing functions, according to the mean value theorem, we 



,Xf„, w „,\A— 1 In, a „,\ *-~ a7T„.A„,A 



(a) 

deduce that J n (v,y) < v x ( 00^^- ) (v Ay). 

First, assume that v + y < 1, then we get J„(w, y) < y/nv x (v V y) A ~ 1 (w Ay) < \/nv x y > 
Next, assume that v + y > 1, then 

Jn(u,2/) < Vnv x {vV y) x ~\v + y) 2X+ %v Ay) <2 2X+e y/^v x {v Ay)(v\/ y) 3 ^ 6 - 1 

< 2 2A+ ^ [(« A y) 2X (v V 2/) 2A+£ l Ae( o,i/2) + (v A y)(v V 2 /) 4A+£ - 1 1a G [i/2,i]] • 
When A e (0, 1/2), we used (v A y) < (v A y) 2X (v V y) 1 ~~ 2X to deduce the bound w A (v A y)(v V 

j / )3A+e-l < /\ y) 2A (-(; V 2/) 1_2A (v V y) 3 ^- 1 < (v A y) 2X {v V y) 2A+,f . 

Thus, in any case 

Jn(v, y) < v^V + 2 2A +^ [(« A y) 2A (« V y) 2A+£ lA e( o,i/2) + (v A y)(v V y) 4A+£ - 1 lA e[ i/2,i]] • 
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Finally, we deduce the bound: 

(2 2A+£ + ^) [(V A yf\v V 2/) 2A+£ l Ae( o,i/2) + {v A y)(v V y)^- 1 !^^] . 
This concludes the proof of Lemma 5.1. □ 
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